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Abstract 



' We introduce a class of isotropic time dependent random fields on the non- homogeneous 

sphere represented by a time-changed spherical Brownian motion of order v £ (0, 1] with 
which some anisotrophies can be captured in Cosmology. This process is a time-changed 
rotational diffusion (TRD) or the stochastic solution to the equation involving the spher- 
ical Laplace operator and a time-fractional derivative of order v. TRD is a diffusion 
i-^H , on the non-homogeneous sphere and therefore, the spherical coordinates given by TRD 

represent the coordinates of a non-homogeneous sphere by means of which an isotopic 
random field is indexed. The time dependent random fields we present in this work is 
therefore realized through composition and can be viewed as isotropic random field on 
randomly varying sphere. 

keyword:random field on the sphere, time-changed rotational Brownian motion, sta- 
ble subordinator, fractional diffusion, Wigner coefficient, Clebsch-Gordan coefficient, CMB 
CN , radiation. 
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1 Introduction and statement of main results 



In recent years a growing literature has been devoted to the study of the random fields on 
the sphere and their statistical analysis. Many researchers have focused on the construc- 
tion and characterization of random field indexed by compact manifolds such as the sphere 
§J = {x 6 M 3 : |x| = r}, see for example [H [231 Ell ES] • In such papers the sphere represents 
a homogeneous surface in which the random field is observed. The interest in studying ran- 
dom fields on the sphere is especially represented by the analysis of the Cosmic Microwave 
Background (CMB) radiation which is currently at the core of physical and cosmological 
research, see for instance |10|ll8j. CMB radiation is thermal radiation filling the observable 
universe almost uniformly |36j and is well explained as radiation associated with an early 
stage in the development of the universe. Prom a mathematical viewpoint the CMB radi- 
ation can be interpreted as a realization of an isotropic, mean-square continuous spherical 
random field for which a spectral representation given by means of spherical harmonics 
holds. Due to the Einstein cosmological principle (on sufficiently large scales, the universe 
looks identical everywhere in space (homogeneity) and appears the same in every direction 
(isotropy)) the CMB radiation is an isotropic image of the early universe [10]. Nevertheless, 
such a nature of the CMB radiation can be affected by anisotropies as those due to the 
gravitational lensing for instance. Recently, a growing attention has been drawn to high- 
frequency or equivalently high-resolution asymptotics for statistics based upon functionals 
of isotropic random fields (see for instance [SJ [23] ) . 

Beside the interest on random fields, particular attention has been also paid, in last 
years, by yet other researchers in studying fractional diffusion equations. These equations 
are related with anomalous diffusions or diffusions in non-homogeneous media, with random 
fractal structures for instance [30]. Starting from the works |16 1 l34 l Rd] much effort has been 
made in order to introduce a rigorous mathematical approach (see for example [33] for a 
short survey on this results). The solutions to fractional diffusion equations are strictly 
related with stable densities. Indeed, the stochastic solutions we are dealing with can be 
realized through time-change that are inverse stable subordinators and therefore we obtain 
time-changed processes. A couple of recent works in this field are \29\ 135] , 

Let D^g(t) with v G (0,1], be the Dzhrbashyan-Caputo fractional derivative of g{t) 
of order v defined by its Laplace transform J °° e~ st D"g(i)dt = s u g(s) — s u ~ 1 g(0) where 
g(s) = j Q e~ st g{t)dt is the Laplace transform of g(t), and D^g{t) becomes the ordinary 
first derivative dg{t)/dt for v = 1. 

Let be a stable subordinator of index v 6 (0, 1) with Laplace transform 



the inverse of the stable subordinator of order v G (0,1). £^ has non- negative, non- 
stationary and non-independent increments (sec [28 J. Let 






We define by 



% = inf{r > : Sf T > t} 






(1.3) 
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be the Mittag-Leffler function. By equation (3.16) in [29] the Laplace transform of 2% is 
given by 

Eexp(-A££) = E v (-Xt v ). (1.4) 
Here and in the sequel we use the fact that i£§J can be represented as 

x = (r sin i? cos ip, r sin $ sin yj, r cos i?) . 

The spherical Laplacian is defined by 

1 



i a 2 i a f . a d 



sin 



G [0,vr], <p G [0,2tt]. 



Sometimes A§2 is called Laplace operator on the sphere. For the sake of simplicity we will 
consider r = 1, that is the sphere of radius one. 

Since the density of Brownian motion *B x (t) on the sphere § 2 started at x G §f solves 
equation 

^jjT^ = d t u{x,t) = A B 2u(x,t), xeSj,t>0 (1.5) 

we say 23 (t), t > is a stochastic solution to (|1.5p . 

Let = cos d(x,y) be the usual inner product in M 3 . With the notation x = 

(sin -& x cos (/Jj; , sin $ x sin , cos $ x ) and y = (sin fly cos </j y , sin fl y sin ^ , cos ■dy ) we have 

cos d(x, y) = cosflx cos $ y + sintfa; sini?^ cos(<p x — ip x ) = (x, y). 

Let Pi, I = 0,1,2, - ■ ■ be the Legendre functions associated to the eigenvalue problem 
A § 2U = —fjiiu. Let {R[, I > 0} be the set of positive real numbers depending on the angular 
power spectrum (defined in section [3|) of no with Rq = 1. The space H s (§1) is a subset of 
£ 2 (§> 2 ) and it is defined in equation (|3.10p . 

Our interest, in this work, is to study the stochastic solution to a time-fractional Cauchy 
problem on § 2 and obtain a new random structure for the sphere by means of which the 
random field is indexed. 

Theorem 1.1. Let v G (0, 1) and s > 7/2. The unique strong solution to the fractional 
Cauchy problem 

d v t u v {x,t\ x ,t ) = A S 2U u (x,t;x ,t ), x, x G §?, < t < t 

2 

u„(x,t ;x ,t ) = u (x - x ), u (- - x ) G # s (§i), 

is given by 

u v (x, t; so, to) = Eno (23*(£^ ) ~ *o) (1-7) 

where *B%,(t), t > to is the time-changed rotational Brownian motion on the sphere 8 2 started 
at x at time to- The explicit solution is written as 

21 + 1 

u u (x, t; so, t ) = Yl ~^T Rl K (* " to ^ Pl ^ x ' ( L8 ) 

l>0 
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This theorem is an extension of the results in [29] and [20] to the spherical Laplacian 
on the sphere. 

%$v°(t) = © X °(-Ct) is called time-changed rotational Brownian motion on Sf. Q3£ (t) is 
a stochastic solution to (jl.6p and has non-independent, non-stationary and non-negative 
increments. Furthermore, for v — > 1, we get that a.s. £^ — >• t which is the elementary 
subordinator (see [7]). When v = 1 the time-changed rotational Brownian motion becomes 
the rotational Brownian motion or Brownian motion on the sphere {*B(t), t > 0. 

We next consider a real-valued random field on the sphere {T(x) : x G Sj}. Let 50(3) 
be the group of rotations in M 3 that can be realized as the space of 3 x 3 real matrices A 
such that A' A = 1$ (where 1% is the three-dimensional identity matrix) and det(A) = 1. 
Suppose that 50(3) acts on Sf with g — > gx. A random field T is said to be n-weakly 
isotropic if K\T(gx)\ n < 00 (n > 2) for every x G Sf and if, for every xi, . . . ,x n G Sf and 
every g G 50(3) we have that 

E[T(xi) x • • • x T(x n )] = E[T( 5 xi) x • • • x T( 5 x n )]. 

We consider the composition of the random field T with an independent time-changed 
rotational Brownian motion QS^(t) G Sf, t > to given by the time dependent random field 

{<F t (x)=T(<S x v {t)),xe§lt>t }, 1/6(0,1]. (1.9) 

We will call {T^(x), t > to} random field on non-homogeneous sphere as the set {58^(t), t > 
to} C Sf is a non-homogeneous set. 

In our view, the composition (|1.9p can be regarded as a random field with randomly 
shifted index. Indeed, we can introduce the shift St such that 

s t x = x + <B° u (t-t ) =»*(*), VxGS? (1.10) 

where = x^r is the North Pole. Thus, the time dependent random field (jl.9p can be 
rewritten as follows 

{^(x)=T( S( x),xG§?,t>t }, */G(0,l]. (1.11) 

Formula (jl.lip says that the compact support of the random field T is affected by the 
random action of the shift (jl.lOp . 

We first obtain in Lemma 14.21 that 

mUaxT = £ J ^ + ^f-^ EKo • • • a in0 ] < 00, G 50(3) 
(is finite) where 

% mj = / T ( x ) Y i* mi (x)Hdx), lj > 0, |m 3 -| < lj, j = 1,2,... ,n 

are the coefficient of T{x) with respect to the orhonormal basis of spherical harmonics {Yi m : 
i,m6Z,l>0, |m| < /} of L 2 (S^). ai jm . are uncorrelated (over /) random variables, A is the 
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Lebesgue measure on the sphere and E[a/ x o • • • a/ n o] is called the angular polyspectrum of 
order re— 1 associated with the field T (See Lemma 4.2 below). The fact that E[T^ {gx)] n < oo 
is a necessary condition for T^(x) to be isotropic. 

The next theorem presents the time- and space- covariance of (jl.9p and this is our second 
main result in this paper. 

Theorem 1.2. For < to < *i < *2 < oo and v G (0, 1] we have that: 

i) forxe Sf, Vg G 50(3), 

27-1-1 

EKW^(P)] = E^ C < E^-ufa-toY) (1.12) 
«>o 

ii) /or i,!/6§j smc/i t/iat x ^ y, 

9/ i 1 

E[S£ (*) S£(j/)] = £ Cj £,(*„, *i, * 2 ) PJ((x, y)) (1.13) 
where C\ is the angular power spectrum defined in equation (|4.4p and 

2 

£i(t ,t 1 ,t 2 ) = Y[E v (-in(ti-to) v )- 

i=i 

A stochastic process t > with E(X(t)) = is said to have short range depen- 

dence if for fixed t > 0, XlhLi E(X(t)X(t + h)) < oo, otherwise it is said to have long-range 
dependence. 

Corollary 1.3 (Time-covariance). For h>0 and v G (0,1], the (equilibrium) time covari- 
ance of the random field %\{x), t > to is given by 

21 + 1 

Cov<z(h;v) :=E[<Z» 0+h (x)<X» t0 (x)} = £ -^C, ^(-/i, fc"), Vg G 50(3). (1.14) 

z>o 

7 -T^Ci> for /i large. (1.15) 

Remark 1.4. From equation (| 1 . 1 5 1) we see that random field %i(x) has long range depen- 
dence for v G (0, 1) since 

oo 

y~] Covx(/i; = oo. 
ft=i 

Corollary 1.5 (Space-covariance). For x,y G 8f suc/i t/zat x 7^ y, t > to > and G (0, 1], 
t/ie space-covariance of the random field %\ (x), t > to is given by 

T v ((x, y),t) := E[S?(x) 3?(y)] = £ ^ C, (^(- W (* - to)")) * fl(fo */»■ (1-16) 
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Remark 1.6. Under suitable choice of the initial condition for the Cauchy problem involving 
the equation ()1.6j) . for all x G §f and t — > to we show that %"{x) — > T(x) and therefore we 
obtain the driving random filed on the homogeneous sphere as expected. Indeed, for t — >■ fo 
we get that T^x) — > T{x) only if the TRD has initial datum uq = 5. Indeed, we have that 

27-1-1 

l>0 
which coincides with (14.61), 



Remark 1.7. Fort — > oo t/ie Mittag-Leffler function goes to zero and thus, formula (|l,16p 
goes to the constant values 

limEffi (x)%Uy)] = ^T 

t— »oo 47T 

where we have used the fact that E u (0) = 1, and Po({%, y)) = 1. For £/ie density of the TRD 
we get that lim t ^ooUu{x,t;xo,to) = 1/4-k and thus, for all v G (0,1), we get a uniformly 
distributed r.v. on the whole sphere. We recall that A(S|) = Air. The covariance of the 
random field %^ x> {x) does not depend on the space and this means that 

Vx G Si, 1^(x) = W(x), Vi/G (0,1) 
where {W(x), x G Sf} is an uniformly distributed noise on the sphere. 
Our third main result is 

Theorem 1.8. For < to < t\ < t 2 < oo and V g G SO(3), we have that 

27-1-1 

EfiKgx^Kgx)} = £ -^-Q exp(- W (* 2 - h)). (1.17) 

l>0 

Theorem 11.81 is not a special case of Theorem ll.2I In this case *B' x (t) = %$f,(t) is a 
rotational Brownian motion on which is Markovian. So we can prove more for this case. 

Remark 1.9. From equation (|1.19p we see that random field l](x) has short range depen- 
dence since 

oo oo , 

Y^mlh^Ux)} = EE^r c < ex P (- W (M) < oo 

h=l h=l l>0 

when Ci ~ l~ s as I oo for 6 > 2. 

By simple conditioning, the independence of and *B l (t), Theorem 11.81 an d Fubini 
theorem we have 

/■OO /"OO 

E[Xi' 1 (flx)Ti' 2 (5x)]= / / E[1 1 Zl (gx)1 1 Z2 (gx)]H(dzi, dz 2 ) 
Jo Jo 

/"OO f'OO 

= J J E[T(^( Zl )T(^(z 2 )]H(dz u dz 2 ) (L18) 
= ^ 21 ^ 1 Cj y y exp(-/j z |zi - 2 |)-ff(^i, dz 2 ) 
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where H{dz\, dz 2 ) = W(£ti-t < z i^t2-t < Using (fTTTKl) and Theorem 3.1 in [21] we 
can deduce the following 

Corollary 1.10. For < to < t± < t 2 < oo and V g 6 SO^), denote by T\ = t\ — to and 

T 2 = t 2 — to we have that 



E^ tl (gx)^ 2 ( 9 x)) = £ 2 ~^Q \e v {-^) + f llT2 — 

i>o n L Jo z 

Remark 1.11. Since by Remark 3.3 in \2Vj we have for fixed T\ = (ti — to) > 
, mi/\ m V f Tl/T2 EJ-fiiT 2 u (l - zY) , 1(1 TV 



dz . 
(1.19) 



jo T%T(1-v)\h T{l + v) 

as T 2 = (t 2 — to) — > 00, we deduce that (x) has long-range dependence. 

Motivations The model we present describe a random motion over a random surface. 
Indeed, if we write the TRD, with staring point (i?o> <Po) a t time to, as follows 

(tft,^) = Q3^°)(i) 5 t>0 

then we have that 

(4^r(^ t ,ft))ei 3 

is a point representing a randomly moving particle over a random surface and therefore a 
motion on a random environment. Apart from the mathematical interest in studying time- 
dependent random fields indexed by a random environment, we want to provide a new 
random field in which the anisotropies of the CMB radiation can be explained. Indeed, 
since the set {58J;(t),t > to} C S| is a non-homogeneous subset of the sphere, in our view, 
the field captures the anisotropies by which the CMB radiation is affected. Furthermore, 
the multiparameter process we present well explain also the observational error due to 
instruments and depending on the observation time. Indeed, the shift representation 
can be considered as a process which well describe such observational error. As time passes, 
we are not observing T{x) but T(x+TRD) where the TRD is a noise depending on time. A 
further remarkable feature is that the new random field T^(x) possesses an angular power 
spectrum, say Ci(t) (depending on time), which goes from an exponential (for v G (0, 1)) 
to a polynomial (for v — > 0) behaviour as the frequency I — > 00. This result plays an 
important role in the asymptotic theory in the high-frequency sense. If the random field is 
Gaussian for instance, then its dependence structure is completely identified by the angular 
correlation function and the angular power spectrum. For non-Gaussian fields, we need 
higher-order correlation functions in order to characterize the dependence structure and, 
in turn higher-order angular power spectra, polyspectra. The asymptotic behavior of (m- 
points) higher-order moments of a random field is therefore of great interest in studying 
the asymptotic nature of the random field, as asymptotic Gaussianity for instance. In our 
case, we do not consider Gaussianity but the high-resolution analysis of reveals different 
asymptotic covariance structure for such a random field, depending on v £ (0, 1) as pointed 
in Remark l4.61 In real data, we get more and more information (or resolution) as / increases. 
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Overview of the paper The plan of the work is as follows. We introduce fundamental 
concepts on fractional calculus in Section [2 In Section [3] we study time-changed diffusions 
on the sphere and establish the connection of the time-changed Brownian motion on the 
sphere to fractional order PDEs involving spherical Laplacian. In Section [U we introduce 
random fields on the sphere and define a new multiparameter process which is the time 
dependent random field on spherical non-homogeneous surface. 

Notations For the reader's convenience we list below some useful symbols: 

• T{x) is the isotropic random field on the (homogeneous) sphere, 

• $)t is the stable subordinator with density h u , 

• £^ is the inverse to a stable subordinator with density l u , 

• B x °{t) is the Brownian motion on the sphere (rotational Brownian motion) starting 
from xq G §1 at to > 0, 

• 23^° (t) = B x °(£t) is the time-changed rotational Brownian motion (time-changed 
rotational diffusion, TRD) with law it„, 

• (x) = T(Q3^(t)) is the time dependent random field on the non-homogeneous sphere. 
Furthermore, 

• or d\ with v G (0, 1) is the Dzhrbashyan-Caputo fractional derivative of order 
1/6(0,1), 

• with v G (0, 1) is the Riemann-Liouville fractional derivative of order v G (0, 1), 
and, for v = 1 

B\ = B] = d/dt. 

2 Inverse stable subordinators and Mittag-Leffler function 

A stable subordinator ft", t > 0, v 6 (0, 1), is (see (7j) a Levy process with non-negative, 
independent and stationary increments with Laplace transform in (jl.ip 

The inverse stable subordinator defined in (|1.2p with density, say l v , satisfies 

Pr{Z u t <x}= Pr{tf x > t}. (2.1) 

Let W£g(t) be Riemann-Liouville fractional derivative of order v G (0, 1) defined by its 
Laplace transform J °° e~ st W t g{t)dt = s u g(s). According to [31 represents a stochas- 
tic solution to 

MB£ + ^ l„(x, t) = 0, x > , t > 0, v G (0, 1) 
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subject to the initial and boundary conditions 

l v (x,0) = S(x), x>0, , 

i,(o,t)=r7r(i-i/), t>o. { > 

Due to the fact that ft", t > has non- negative increments, that is non-decreasing paths, 
we have that £^ is a hitting time. Furthermore, for v — > 1 we get that 

lim S) y t = t= lim £% 

v — y 1 V— >1 

almost surely ([7]) and therefore i is the elementary subordinator. From the relation (|2.ip 
and the Laplace transform (jl.ip . after some algebra, we arrive at the Laplace transform of 
£ u given by (H3D- 

Next we state some of the properties of the Mittag-Leffler function. Let v £ (0, 1] As 
we can immediately check E v (0) = 1 and (see for example [151 EZ] ) 

< E v (-z v ) < —— < 1, z £ [0, +oo). (2.3) 

L ~\~ Z 

Indeed, we have that 



whereas, 



Thus the Mittag-Leffler function is a stretched exponential with heavy tails. Furthermore, 
we have that (see [151 formula 2.2.53]) 

W z E u (fiz u ) = + iiE v (iiz v ), fi £ C. (2.6) 

Dzhrbashyan-Caputo derivative of order v £ (0, 1) is defined in [9] by 

D ^=f(Tb)r^ (z - sr " <is - ( "» 

Formula ()2.7p can be also written in terms of the Riemann-Liouville derivative as follows 
(see [USEE]) 

££/(*) = lD£/(s) - /(z)|, =0+f ^-^, (2.8) 



and therefore, formula ()2.6[) takes the form 

D v z E v (iiz v )= fiE^z"), neC, i/G(0,l). (2-9) 

Hence in this case we say that E u (^iz u ) is the eigenfunction of the Dzhrbashyan-Caputo 
derivative operator D v with the corresponding eigenvalue fi £ C. 
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3 Time-changed Rotational Brownian motion on the sphere 

In this section we define a measurable map from the probability space P) to the 

measurable space (§ 2 , i3(§ 2 ), A*b) which is the time-changed rotational Brownian motion 
t > to, v G (0, 1] with starting point xq G Sf at time to > 0. Such a process can be 
regarded as a time-changed rotational diffusion (TRD) or a rotational Brownian motion on 
the sphere § 2 time-changed by an inverse to a stable subordinator. Thus, the composition 
we deal with is written as 53^° (t) = B x °(£,^), t > to where B x ° is a rotational Brownian 
motion starting from xo G §f at to > and the time-change is given by t > which is 
the inverse to a stable subordinator of index v G (0, 1). The TRD represents the spherical 
counterpart of the fractional diffusion on bounded domain driven by the fractional equation 

{D\ - A)u(x, t) = 0, x G W 1 , t > 
n(x, 0) = uq(x) 

where A = Y17=i * s ^ ne Laplace operator and whose solutions can be written in terms of 
Wright or Fox functions ( I13|, [35] ) and by using time-changed Brownian motion in M. n ; see 
for example |28t |2"5] . In [55] the authors presented the explicit one dimensional solutions for 
some particular order a of the fractional derivative. For information on fractional diffusions 
on unbounded domains the reader can consult, for example, the works [XT], [221 HH 155] . 
In [3] the fractional Cauchy problem involving an infinitely divisible generator on a finite 
dimensional space has been investigated while [29] studied the fractional Cauchy problem 
(|3.ip in a bounded domain D C M™. In [29], the authors found that (in our notation) the 
strong solution of (|3.1|) in a bounded domain D with Dirichlet boundary conditions is given 

by 

u(x,f) = Euo(B x (£H) 1 (Tb(B )>£ ? ), x G V, t > 

where 

r„(B) = inf{s > : B(a) £ P} 

is the first exit time of Brownian motion B from T> C M n . 

Our aim, in this section, is to study the solution to the fractional Cauchy problem 
involving the spherical Laplacian operator which extends the results in [20} ?] and prove 
Theorem 11.11 

Remark 3.1. We say that A § 2U exists in the strong sense if it exists pointwise and is 
continuous in D. 

Similarly, we say that f(t) exists in the strong sense if it exists pointwise and is 
continuous for t G [0, oo). One sufficient condition for this is the fact f is a C 1 function on 
[0, oo) with \f'(t)\ < ct 7_1 for some 7 > 0. Then by (|2.7|) . the Caputo fractional derivative 
d^fit) of f exists for every t > and the derivative is continuous int>0. 

It is well-known that ([H]) the solutions to the eigenvalue problem 

Ag2ti(x) = —/j,u(x), x G S 2 (3-2) 
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is solved with a sequence of eigenvalues fii = 1(1 + 1), I £ Z, I > with their corresponding 
eigenfunctions given by the spherical harmonics (here we use the fact that x € can be 
represented = (sin t? cos (p, sin $ sin cos i?) ) 



21 4-1 (I — 77? V 

Hm(aO = *lm(0,¥>) = \/ 4 ^ a + 7n)! P ' m(c ° S ^ )e ^' H " *' * G t '^' ^ G [ °' 27t] 

(3.3) 

(or linear combination of them) where 

p lm ( Z ) = (-ira- f)™/*—F\{z) 

are the associated Legendre functions, and the well-known Legendre polynomials Pi are 
defined by the Rodrigues' formula 

Remark 3.2. It is important to note here that the operator considered in Theorem 3.2 in 
f2U\l for the Jacobi case (the case a = b = in their notation)is different from equation 
(|1.6p and ()3.2[) here. First, their fractional Pearson diffusion equation is in K whereas our 
equation is in . Second, the spherical harmonics in (|3.3|) are different from the eigenfunc- 
tions called Jacobi polynomials in their notation. Jacobi polynomials are the eigenfunctions 
corresponding to the eigenvalue problem 

n 2^d 2 9{x) dg(x) 

(1 — x ) — —5 2x — - — = —Xg(x),x e R 

dx z dx 

and it has eigenfunctions gi(x) = (— V) l 2 l V.Pi(x) and corresponding eigenvalues Xi = 1(1 + 1). 

We list some properties which will turn out to be useful further in the text: for all x£§J 
(symmetry) we have that 

Y? m (x) = (-lriWx), (3-4) 

and in spherical coordinates <f>) = ^Im($> ~t°) where * stands for the complex conju- 

gation; for all l\,l-2i mi,m2 (orthonormality) we also write the formula (|A.9j) as 

YhrnAx)Y i ; m2 (x)\(dx) = 6 l ™ (3.5) 

I si 

where S l £™ 2 = 5 l £ 5™^ are the Kronecker's delta symbols (|A.8p and A(-) is the Lebesgue 
measure on E>f given in spherical coordinates by 

\(dx) = s'm'dd'ddip for x = (sin #cos (p, sin -i?sin ip, cos $); 

for all x, y G (addition formula) we have that (see page 339 in p3] for this and other 
properties of spherical harmonics) 

+i 9/4-1 

Y, YUx) YC m (y) = ~^-m^ V)) (3-6) 
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where (x,y) = cos d(x,y) and d(x,y) is the usual spherical distance on Sf. Furthermore, 

£ y, m (x)l^(x) = ^ (3.7) 

m=— I 

and Pi((x,x)) = 1; for all cc,y G §f (reproducing kernel) the following holds 

47T 



P il ((x,z))P i2 ((z,y))A(dz) =^- I __p /i ((x,y)). (3.8) 



/§2 

Let / G £ 2 (§f). Since {Yj m : Z G N,m G Z and |m| < Z} is dense in L 2 (§f). The function 
/ (see also the Peter- Weyl theorem on the sphere in [27] and the references therein) can be 
written as 

oo +1 

f{x) = f(0,<p) = ^ E *lm{f)Yl m {p,<p) 
1=0 m=-l 

which holds in the L 2 sense, where 



aim(f)= m<p)Y? m (0,<p)Bm0d#dtp. 

We define 

+i 

Mf) = E i a ^(/)i 2 ( 3 - 9 ) 

m=—l 

as the angular power spectrum of /. The angular power spectrum can be considered 
as a power law of index, say w, which can be different from zero (red spectrum for w < 
and blue spectrum for w > 0) or not (white spectrum). If Ai(f) decays exponentially as 
Z — > oo, then / is a real and analytic function and turns out to be infinitely differentiable 
on the sphere. In general, for s > 0, we introduce the Sobolev spaces 

H s ' 2 (Sj) = j/ G L 2 (S 2 ) : ||/|| Sj2 = fj(2Z + l) 2s Mf) < oo J . (3.10) 

The Sobolev space ()3.10p is the closure of the set of all spherical polynomials with respect 
to the norm || • || s 2. Furthermore, H S,2 (S 2 ) is a Hilbert space, with inner product 

00 +z 

(fl,f2)s,2 = Y,( 2l + 1 ) 2S E a lm(h)aL(f2), fl, h ^ H**{§\). 
1=0 m=-l 

Obviously, #°' 2 (§?) C L 2 (Sj) and, for s' > s, we get that 

/ g # S '> 2 (S 2 ) => /GiF> 2 (S 2 ). (3.11) 

We write H S {§1) instead of H s ' 2 (Sj) and say that / G H s (Sj) is infinitely differentiable 
provided that A\ ~ l~ 2s ~~ 8 with 6 > 1 as Z — > 00 (s > 0) by the Sobolev embedding 
Theorem. 

We are now ready to give the 
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Proof of Theorem 11.11 The proof follows the main steps in the proof of Theorem 3.2 in 
|20| and Theorem 3.1 in [29]. Since operator ASf is quite different operator than those in 
these two papers, we give all the details of the proofs. The proof is based on the method 
of separation of variables. Let u(t,x) = G(t)F(x) be a solution of (|1 .6j) . Then substituting 
into (|1.6p . we get 

F(x)d^G(t) = G(t)A sl F(x). 
Divide both sides by G(t)F(x) to obtain 

8?G(t) _ A S 2FQr) 
G(t) ~ F(x) ~ * 

Then we have 

d%G(t) = -/xG(t), t > (3.12) 

and 

A § 2F(x) = -aiF(x), x e Sf. (3.13) 

By the discussion above, the eigenvalue problem (|3.13p is solved by an infinite sequence of 
pairs {nim = 1(1 + 1), Y/ m : I > 0, — I < m < +/} where the spherical harmonics Y[ m forms 
a complete orthonormal set in L 2 (§>f). In particular, the initial function no regarded as an 
element of L 2 (§> 2 ) can be represented as 

oo I 

= E E b i m (xo)Y lm (x). (3.14) 

1=0 m=-l 

where bi m (xo) = J s2 uo(x — xo)Yj ¥ m (x)\(dx) = RiY* m {xo). By equation (|2.9|) we see that 
aim(t;xo,t ) ■= bi m (xo)E u (-m(t - t ) u ) solves f|3.12j) . Sum these solutions bi m E u (-fn(t - 
to) u )Yi m (x) to dL6]), to get 

oo I oo I 



u u (x,t;x ,t ) = bi m (x )E u (-fj,i(t - t ) u )Y lm (x) = ^ ^ a im (t; x , t )Yi m (c 

1=0 m=—l 1=0 m=—l 

oo +1 

= E E ^ Y ^°) E v {-Hi{t- to)") *W*)- 



Z=0 m=-Z 



(3.15) 



By the addition formula (|3.6p . we can write 

2/ + 1 

-u„(x,t;2;o,to) = ^2 Ri E u (-fn (t - t ) u ) — Pi((x,x )). (3.16) 

This is similar to the form of the solution of a Factional Jacobi diffusion that was worked 
in section 7.4 and 7.5 in 131 1 for the case a = b = in their notation. 



13 



It remains to show that (|3. 15[) solves (|1.6p and satisfies the conditions of Theorem 11.11 
Since no € H s (§1) we know that 

l l 

Mu ) = Yl (blm(xo)) 2 = £ ( R l Y *m(^)) 2 < C(2l + l)~ 2s (3.17) 

m=—l m=—l 

where C = ||uo|| s ,2- The angular power spectrum for u u (x, t; xo, to) in (|3.15p is then given 
by 

Ai(t,t )= £ \ai m (x;xo,t )\ 2 

\m\<l 



\m\<l 



(3.18) 



E^-m (t - t ) u )f Mu ) 
<Mu ). 

Hence for fixed t > to a s I — > oo, using (|2.5|) we can see that 

Mt,t ) ~ (ju(t-to) v )- 2 M*o) ~ ( (t-to)^!-^ ) r4 ^(«o)- ( 3 - 19 ) 

Step 1. First we check that 

||tv(-,t;a;o,*o)|||2( S 2) < oo 

uniformly in t > to- First we notice that From the Parseval's identity, Equation (|3.18p and 
Equation (|2.3|) we get that 

\\u v (;t;xo,to)\\ 2 L2( p 2 ) =E \ a im(t;x ,t )\ 2 = y~]Ai(t, t ) 

lm l>0 

<^2M u o) = IM- - x o)\\h(§D 

l>0 

this proves that u v € L 2 (S 2 ). 

Step 2. Next we show that the series (|3.15|) is absolutely and uniformly convergent for 
t > to > 0. 
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Using (|2.3p . (|3.17p . (|3.6p and Cauchy- Schwartz inequality we have 



oo +1 

Y E \ a lm(t;X ,t )\\Yi m (x)\ 
1=0 m=-l 
oo +£ 

< J]£ \bim(xo)\E v (fxi(t-t Q ) w )\Yi m (x)\ 

1=0 m=-l 



oo 
i=0 



+/ 



E (^(w(*-*o) I ')) 2 |^(^o)P 



m=—l 



m=—l 



< (t-t )-^r 2 Vc(2/ + i)- s A /^^< 



i=0 



47T 



(3.20) 



since s > we have s + 2 — 1/2 = s + 3/2 > 1. Hence the function in (|3.15p is absolutely 
and uniformly convergent. 

Step 3. We next show that d!^u u (x, t; xq, to) exists pointwise as a continuous function. 
Using pH Equation (17)] 



dE v {-in{t-toY) 



dt 



<c Mt -, tor :l. <cMt-t y-\ 



i + ^{t-toY 



and (pTTTp and t^Ffl) we get 

oo 

\d t u v (x,t;x ,t )\ < X \ b lm(x )\ 



n=l 



dE v {-m{t-to) v ) 



dt 



\Yl m (x) 



oo +1 
1=0 m=-l 

:= c^-to)^ 1 ^) 



x 



< dt-toy-'Y 



/=0 



+/ 



E i 6 '-( 



m=—l 
oo 



m=— i 



< cy/UJ^(t - t ) y_1 E^ + X ))( 2/ + 1)" S+1/2 < oo 



z=o 

since — s + 5/2 < —1 when s > 7/2 = 1 + 5/2. Hence bv Remark 13. 11 d^u u (x, t; xo, to) exists 
pointwise as a continuous function and is defined as a classical function. 

Step 4. We show here that we can apply A S 2 and term-by-term to the series (|3,15p . 
For this we need to show that the series (|3.2ip is absolutely and uniformly convergent for 
t > t > 0. 
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From (|3.2p . term by term we get 

oo +i oo +1 

Y Y a lm (t;x ,t )Yi m (x) = - ^ Y a im(t;xo,to)fiiY lm (3 



1=0 m=-l 



1=0 m=-l 



11 + 1 

W ^ E u {-m (t - t ) v ) — Z-i^fo x )) 



(3.21) 



l>0 



4tt 



By the previous step, using (|2.3|) . ()3.7[) and the Cauchy- Schwartz inequality we get 

oo +1 

i=0 m=~Z 
oo +Z 

< J2J2 \ b i m {xo)\E v {^{t-t Q Y)ixi\Y lm {x)\ 

1=0 m=-l 



< Y^\E u {^{t-t Q Y)^ 



< 



1=0 

oo 

El 

1=0 



+1 

Yl \ b lm(xo)\' 

m= — l 



m=—l 



C(2l + 1) 



(21 + 1) 

47T 



OC 



(3.22) 
(3.23) 



l=o 



since s > 3/2 we have s — 1/2 > 1. Hence 

oo +1 oo +Z 

Y Y 9 t a im(t;x ,t )Yi m (x) = Y Y a im(t;x ,t )A § 2Yi m (x) 

1=0 m=—l 1=0 m=—l 

since the two series are equal term-by-term, and since the series on the right converges 
absolutely. 

Now it is easy to check that the fractional time derivative and A§2 can be applied term 
by term in (|3. 15|> to give 

(dt ~ A § 2)u u (x,t;x 0: t ) 

oo +1 

= Yl Yl hm{xo)[Yi m {xWt E u {-tn{t - tof) - E u {-^{t - t Y)A §2 Y lm (x)] = 0, 

1=0 m=-l 

so that the PDE in (|1.6p is satisfied. Thus, we conclude that u defined by (|3.15p is a classical 
(strong) solution to (|1.6p . 

Step 5. We show that *B x (£(t)) is a stochastic solution to (|1.6p . The solution to dtu\ = 
A§2iti subject to the initial condition uq is written as ui(x, t; xq, to) = Euo (Q3 X — to) — %o) > % G 
Sf, t > t > or, in terms of the convolution semigroup, as 



ui(x,t;x ,t ) = P t -t u (x - xq) 
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where the strongly continuous semigroup Pt on L (Sf) is satisfies 

PtUl-Ux 

lim = ZX532U1 

Since A§2 is self-adjoint operator by Remark 3.34 in [27] . We define the operator 



P?f(x)= / l u (s,t)P s f(x)ds. 
Jo 

Finally, we obtain the stochastic representation of the solution. Since {Y/ m } forms a 
complete orthonormal basis for L 2 (S 2 ), the semigroup 

00 1 

E[u (<B x {t - t ) - so)] = P t -t u (x -x ) = Y, Yl e-^b lm {x )Y lm {x) 

l=Q m=-l 

is the unique solution to (|1.5p where 

P t -to uo(x - xo)Y l * m (x)\{dx) = e'^kmixo) 



we can write by using Fubini Theorem together with equations (|1.4|) and ()3.15p to get 

00 1 00 1 



1=0 m=-l 



u u (x,t;x ,t ) = ^2 E h im{xo)E v {-Hi{t - t Q ) v )Yi m (x) = ^2 Y a im(t;xo,to)Y lm (a 

1=0 m=-l 

00 j /*oo 

= E E b i m (x )Y lm (x) J e-^H v (t,y)dy 



1=0 m=-l 



00 I 



lu(t,y)dy 



^ k m (x )Yi m (x)e 

,1=0 m=-l 

i 

P y -t u {x - x )lv(t, y)dy = P"_ tQ uo(x - x ) 

E[uo(<S x (y)-x )]l v (t,y)dy 
= Eiuo^W-J)}- 

Uniqueness follows by considering two solutions u±,U2 with the same initial data, and 
showing that u\ — 112 = 0. 

□ 

Remark 3.3. We notice that if the initial condition in (|1.6|) is the Dirac delta function 

u (x - xq) = 5(x - xo) 
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then, from the completeness relationship for spherical harmonics 

n +1 

5(x-x ) = lim 5 n (x-x ) = lim V] } Y* (x )Yi m (x) 

Z=0 m=-Z 

we can lonie i/iat 

aj m (*;x ,to) = Y l * m (x )E v (-fj,i (t - to) u ). 
As we can check, this implies that Ri = 1 for all I > wni/i 

2/ + 1 

Ai(t,t ) = — , for all t,t >0 

47T 

and also that 

21 + 1 



<UO) = x; 



47T 
Z=0 



Remark 3.4. We check that the transition density (|1.8p integrates to unity. We recall that 
Yqq(x) = l/y/4TV for all x G §f . Furthermore, from the integral ()A.9p , we /lave i/iat 



y Zm (x)A(dx) = / ^ m (x)y o(^)A(dx) = 

From this and the fact that E u (0) = 1, we obtain 

u v (x,t;xo,t )X(dx) = / u u (x,t;x ,t )\(dx ) = R 



which, under the hypothesis of Theorem \l.l[ leads to the claim. 

Remark 3.5. For xq = xn, that is the North Pole, we have that (see Section 5.13.2 of 
ffl) 

Yim{xo) =0 form^O and Y l0 (x ) = i/ ^^ 1 
whereas (JJDj), for x = tp) G Sf , 



A /^!p,(costf). (3.24) 



T/ie distribution (11,81) becomes 



97 + 1 

*; rr^v, to) = X] I - ^ " *o)") fl(cos 1?) (3.25) 



47T 

Z>0 



/or < 1? < vr and i > t > 0. 
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Remark 3.6. For v — > 1 we have that (see J^j) £j -4-' t and therefore the inverse process 
£}, t > 0, becomes the elementary subordinator t, t > 0. On t/ie ot/ier hand the fractional 
equation of the Cauchy problem (|1 .6j) . /or v = 1, ta/ces t/ie following form 



d t ui = A S 2Ui 

mtt solution, from the classical Sturm- Liouville theory, given by 

+i 

ui(x,t;x ,t ) = ^2 ^2 Rl ex P _ *o)) Ylm(x)Y* m (x ) (3.26) 

;>o m=-l 

which is in accord with the property of the Mittag-Leffter function E\{—z) = e~ z . Formula 
(|3.26p can be also written as 

ui(x,t;x ,t ) = Eiio(23 x (t - to) - %o), x,x G Sf, t > i 

where *B x (t), t > is the rotational Brownian motion or Brownian motion on Sf started at 
x at t = to- 

Remark 3.7. For ^ ^ 1 and no 7^ 5 the following holds true 

u u (x 2 ,t 2 ;x ,t ) / / u I/ (x 2 ,t2;xi,ti)n ! ,(xi,ti; sco,io) X(dxi) (3.27) 

where to < ii < t 2 and xq,xi,x 2 £ Sf. Indeed, from (|3.8|) we /tawe t/jat 

; 47T 

^ 1 ((x 2 ,xi))Pi 2 ((xi,x )) A(dxi) = (5^ — - - /^((a^so)) 



and therefore 



u u (x 2 ,t 2 ;x 1 ,t 1 ) u u (x±, ti; x , to) A(drci) 



= Z) ^j^ 12 * ^(-W(*2 - tiD Evi-mfa - to)") Pi((x 2 ,x }). 

l>0 

which coincides with u u (x 2 ,t 2 ;xo,to) only if is = 1 and no = 5. Indeed, for v = 1, 

E v {-^i{t 2 - t x ) v ) E v (-m(ti - t n = E v {-^{t 2 - t y) 

and, for uo = 5, Ri = 1 /or aW £ > 0. 
Remark 3.8. We observe that 

P{d(x, aS*(t)) < p} = P {rJ (a . p) (»S) > t > to} 

w/iere 
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and 

B(x,p) = {y G §> 2 r : d(x,y) < p}, p€(0,ir). 
From the fact that <BJ = 53^ for all g G SO(3), we can also write 

or equivalently that 

T/ie above distribution (see for example J2J/J solves the Dirichlet problem 

A§2 q{&, ip; 0„, Vo) = 0, < < tf < 7T, V?, 93 G [0, 2vr) 

g(t?o,^;^o,^o) = %> - Vo), 

and is written as 



(3.28) 



1 COS 1?n — COS 1? 

O , <A) = ^ ; ; ° : , 7 r- 3.29 

Z7T 1 — cost/ cos^o — sin 77 sm Vq cos(9? — <po) 

The solution, say u u , to the problem (jl.6p on the bounded domain 

V = {(d,<p) : < $ < $o < ir, (p,tpo £ [0.27T), OWo) 

can &e written as 

u u {x,x ) = En (<S x m)) l( r o ((S o )>t _ 4o ), x G Z>, t > t ? (3.30) 
5ee fUP) / /or more on fractional Cauchy problems in bounded domains. 

4 Random fields indexed by spherical non-homogeneous sur- 
faces 

4.1 Random fields indexed by the sphere 

In this section we will consider a n- weakly isotropic random field {T(x); x G Sf} on the 
sphere for which 

KT(gx) = 0, and T{gx) = T(x), for all g G 50(3) 

here we recall that = means equality in distribution of stochastic processes and SO (3) the 
special group of rotations in M 3 . 

As already pointed out in Section [31 the triangular array 

{Y lm : I > 0, m = 
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of spherical harmonics represents a set of eigenfunctions for the Laplacian on the sphere 
and an orthonormal basis for L 2 (§ 2 ). Hence we can see that 2- weakly isotropic random 
fields admit the spectral representation 



l>0 l>0 \m\<l 



where 



a h 



[ T(x)Y l * m (x)X(dx) 



(4.1) 



(4.2) 



is a set of random spherical harmonics coefficients that are Fourier random coefficients and 
A is the Lebesgue measure such that A(§ 2 ) = 4-7T. In formula (|4.2p the symbols Y^ m stands 
for the complex conjugation of Yj m and convergence in (|4.ip holds in the mean square sense, 
both with respect to L 2 (dP) for fixed x and with respect to L 2 (dP <S> X(dx)), i.e. 



and 



lim E 

L— »oo 



( L +l X 

lim E T(x) - y~] a lm Yi m {x) 

L— >oo \ *■ — ' ^ — ' / 

\ «=0 m=-l / 



L +1 \ 2 

T{x) - ^2 ^2 ai m Yi m (x) \{dx) 

1=0 m=-l ) 




0. 



Under isotropy, {a/ m } are zero- mean valued and uncorrelated over I and m with E \ ai m \ = 
C\ where C\ is the angular power spectrum of the random field T. We notice that (from 
the isotropy of T) C\ depends only on I as shown in [61 [26] . Furthermore, for the harmonic 
coefficients associated with the frequency I, we have that (see formula (|3.4p ) 



a lm = (-l) m ai 
and the following orthogonality property holds 



^[a hmi ai 



2m 2 l 



(5^ is the Kronecker symbol (|A.8|) ). From (|4.4p we immediately get that E|t% 
I > 0. As in ()3.10p . for the angular power spectrum C\ we have that 

^2^2 (21 + l) 2s E|a Zm | 2 = ^2(21 + l) 2s (2Z + 1)Q < oo 

l>0 \m\<l l>0 



(4.3) 

(4.4) 

= Q, 

(4.5) 



for some s > 0. In the light of the previous discussion, for the zero-mean n-weakly isotropic 
random field T we can write the following expression for the covariance function 



07 I 1 

E[T(x)T(y)]=J2-[ —QPidx^y)), x,y€Sl 



l>0 



4tt 



(4.6) 
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where (x, y) is the angle between x and y. Indeed, from the representation (|4,ip . we have 
that 

E[T(x)T(y)] = ^ E E i T h ( x ) T h (V)] ( 4 - 7 ) 

h>0l 2 >0 

where 

+h 

T h( z )= E a hmY hm (z), i = 1,2 

m=— li 

and, by taking into account formulae (|4.3p and (|4.4p . we get that 

E[T ll (x)T l2 (y)}= ]T £ ft C ^i mi WC 2 (2/) 

|mi|<Zi |m2|<Z2 

=C h £ Yhm^yilmM- 
\mi\<h 

From the addition formula (|3.6p the covariance (|4.6p immediately follows. 

By considering that Pi({x,x)) = 1, from (j4.6|) and (j4.5|) we also obtain that 

E[T(*)] 2 =£^-Q<oo. (4.8) 
z>o 

Under isotropy, the harmonic coefficients a\ m are zero-mean and uncorrelated over I. Fur- 
thermore, if T is Gaussian, then a/ m are Gaussian and independent random coefficients (see 
[B]). We will use throughout also the following fact (sec [26J) 

EK mi ...a lnm J = (-ir Y C ^Xi^ P ^n(^---^n-^ (4.9) 

Ai...A n _3 

where C, Al '," A ™r 2Z ™ TV , m ™ is a convolution of Clebsch-Gordan coefficients and Pi. i (■) is the 

iimi...l„-im„_i ii,—,in\ J 

reduced polyspectrum associated with the isotropic random field T. If n = 3 for instance, 
then the formula (|4.9|) becomes 

E[a hmi a l2m2 a l3m3 ] = {-l) ms C l ^ m P lMa . (4-10) 

For more information on random fields on the sphere the reader can consult the book 
by [27] whereas, we refer to the book by [1] for random fields and geometry. 

4.2 Random fields indexed by the non-homogeneous sphere 

In this section we focus on the main object of the work which is the time-changed (we mean 
composition of processes) random field 

T t (x)=T(*8 x v (t))=T{x + <B(£U )), t>t >0, xeSl ue (0,1] (4.11) 

where is the time-changed rotational Brownian motion starting from x at time 

to and {T(x), x G Sf} is the (Gaussian) random field depicted in the previous section. We 
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assume that T is independent from 23*. Random field with randomly varying parameters 
have been also studied in the interesting work [2]. Based on the representation (|4.ip of T, 
for the random field 

{% v t {x) : x £ S?, t > t > 0}, ue (0, 1] 



/ 



we write 

c 

St (*) = E E a ta >W23*(i)), t > to, x &§i, v £ (0, 1] 

i=0 m=-Z 

which must be understood in the L2 sense, that is 

L +1 



(4.12) 



lim E 

L— >oo 



E 



-EE a lm Y lm (K(t)) 



1=0 m=-l 



- 2 









0. 



We recall that the TRD 23*(t), t > t is a measurable map from (fi,&8,P) to (S?, #(§?), A®). 
In the rest of the paper use the double summation notation whenever appropriate 

00 1 

£ = ££• 

lm 1=0 m=—l 

Remark 4.1. Let us consider the random shift St ■ x 1— )• 03^ {€), such that 

s t x = x + r&l{t-to)=re>l{t) 
where = xjv 2s the North Pole. The process (|4.1ip can 6e therefore regarded as 

T^x) = T(six), x £ E>1, t > to 



(4.13) 



(4.14) 



where the random field T is indexed by stx and therefore under the action of st- Formula 
(|4.13p says that the spherical coordinate x £ §>f is affected by some noise depending on time. 
As we will see further in the text, for v — > 1 we get that 23^ — > 93° and therefore the noise 
becomes the standard rotational Brownian motion. For v £ (0, 1), the random environment 
in which T is indexed is represented by the TRD 23* and therefore is a non-homogeneous 
sphere. 

As expected, the properties of (|4,lip are strictly related with those of T. In particu- 
lar, we begin our analysis by presenting the following result concerning the higher-order 
moments of the composition %^{x). 



Lemma 4.2. Let the previous setting prevail. The following holds 

E[^{gx)] n = E[T(x)] n , t>t >0, v £ (0, 1], re £ N g £ 50(3). 



(4.15) 
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E{[£a lro 15m(®S(*)) 

^ lm 



Proof. The claimed result follows by observing that 
E[%\x)] n =E 

=E 
=E 

From the fact that (see |26j) 



^{[Y^ai m Y lm {y)\ n \y = ^ x At) 

. ^ lm 

E{[T(y)Ay = <B*(t)}], Vx G Sf. 



Vx G S? 



E[T(y)] = W ^— E[a il0 ---a in0 ] 

which does not depend on y G as formula (j4.9[) entails, we obtain that 

E[S?(aO] n =K[T(y)]", Vy G S? 

which, in turn, does not depend on x G Sf. Furthermore, 03^ = for all g G 50(3) and 
this concludes the proof. □ 

Proposition 4.3. For x G Sf , t > to > and G (0, 1], we have that 

EY lm (K(t)) =Eu(-Hi(t ~ to) v )Y lm (x) (4.16) 

and 

EYMXXt)) =E v { rl i l {t - t Y)YC m {x). (4.17) 



Proof. We recall from Remark 13.31 that, for v G (0, 1], 



+/ 



,(y,t;x,t ) = J2 E A-M(-t-to)l Yl Y i™(y)Yi m {x) 



1=0 



m=—l 



(4.18) 



y, x G §i, t > to > 0, is the probability density of the TRD 5S^(t), t > to- From the property 
(|3.4p and the orthogonality of spherical harmonics ()3.5|) . we get that 



F^(y)H ro (y)A(dy)=^ 
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and thus 

E[r, m (»£(t))] = f u u (y,t-x,t Q )Y lm (y)\(dy) 

oo 7 



2 2 ^(-^(t-tor)l^(a;)(-ir/ 2 ^(y)^ m (l/)A(dy) 

7 =0ft=-7 



7=U ft=— 7 

oo 7 



=E E ^(-^(*-*o) !/ )(-ir^(^)^^ 

7=0 ft=— 7 

= [by (1331)] = E„(- w (t - to)")*^)- (4.19) 
Formula ()4.17|) can be obtained from result (|4.19|) by considering the fact that 

E[^(2£°(t))] = (E[r lm (»*°(t))])* = (« - ioD^o) (4-20) 

and the proof is completed. □ 
Proposition 4.4. For x G Sf , t > to > and z> G (0, 1], we Ziaue ZyW 

E [^(BSCt)) y^(!8J(t))] = E C 7% ^H*r(* " W^OO ( 4 - 21 ) 

7,K 



. 7jK -(-l) (2^ + 1)^/^—1 Q Q ) 1 k 772 — Tfl 



is written in terms of Wigner 3 j -symbols (see Appendix^). 
Proof. From (|4.18p we write 

E[Y lm (K(t))Yr m (K(t))] 

u v (y, t; x, t ) Yi m {y) Y* m {y) X(dy) 



•2 

oo 7 



E E K(-^(t-t r)Y; K (x) Y 1K (y)Y lm {y)Yf m (y)\{dy) 

7 =0k=-7 

oo 7 . 

(-^E E ^(-^(i-io)")^^) Y 7K (y)Y lm {y)Yi- m (y)X{dy). 



7=0 K=— 7 

From the fact that (see formula (jA.10|) ) 

Y JK (y)Y lm (y)Y l „ m (y)\(y)=¥ 

2 ' 



7 Z I \ ( 7 Z Z 



7 V 0/\k 771 — m 



where 



^ = (21 + 1)^^, (4.22) 
we arrive at formula (|4.2ip . □ 
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Remark 4.5. We use an alternative approach in order to show that 

Var[% v t {x)) = Var[T(x)} = ^ ^j-^C, < oo 



l>0 



or equivalently that 

E[T(s t gx)] 2 = E[T(x)} 2 , V<? G 50(3), t > t > 
as already stated in Lemma \4-S\ 
Proof. We have that 

Im 



E[%\x)Y =E 



=E 



where 



[E J2 ^ mi a l2m2 Y hmi (^l{t))Y l2m2 {^l(t^ 
hmi hm2 

E[a, imi o laroa ] = (-lr^-^a 



5<8 



and (-l) m l^ m = K* Thus, we obtain 



E[1»(x)] 2 =E 



iimi 



where 



and 



C 7 ,„-(-l) ( 2Z + 1 )V — 4^ — V j 1 , m -m 
as shown in Proposition 14.41 From this, we write 
+/ 



•i=—i 
+i 

£ (-ir j> 7 Eu-Mt-tdrw 



m=—l 



7,k 



7 Z Z \ / 7 Z Z 



/ \ k m — m 



£ ^ ^(-^r(t - OTW ( o ) B" 1 )™ ( I ' ' 



m — m 



(4.23) 



: £ C Zl E [^ imi (^(t))F z * mi (^(t))] 

iimi 



E^^^^^d^C^ ^(-^(i-OTW (4.24) 

7K 
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where 



From (fOfl) and the fact that (00 



1 



we get that 



7 I I 




I I 
, 

+1 



) E (-^ ( 

J m=-l 



7 I I 
K m —m 



5°5° 



and, by taking into account the fact that (see for example |40j ) 

1 



*oo(tf,0) 



2^ 



and E u (0) = 1, i?o = 1, we obtain that 



I III 



7K, 



^ (x)^C l ^ = E[r(.)] 2 , 



4tt 



which is the claimed result. 



□ 



In order to study the time- and space-covariance of the random field 

{5£(x) = T(®S(t)); x G S?, t > t > 0}, i/ G (0, 1]. 

We assume that for given starting points x,y G Sf at to > 0, the processes QS^(ti) and 
33^(t2) are independent if and only if x 7^ y. If x = y, then SB^(ii) and 23^(i2) are two 
dependent copies of a time-changed rotational Brownian motion starting from x = y at to. 
From the representation 

00 +1 

%{x) = T(<8 x v (t)) = E E * > *o (4.25) 



we can write 



EpZWSUv)] =E 



;=0 m=-2 



e^ [r(^(t 1 ))r(^(t 2 )) 



#<8 ® 5<B 



=E 



EQ^(^(tl))^(^(t 2 )) 



Irn 



&8 ®#<8 



(4.26) 
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for some t\, t 2 > and x,y G Sf . Next we give 

Proof of Theorem 11.21 We first consider the starting point y = x at time to > 0. 
i) Due to the fact that 

P{®*(t 1 ) G dy,<B x M G dx} =P{S^(to) G dx}P{»J(t 1 ) G dj/|«8^(to) G cfe} 

-»{®S(tl)€dy|»*(to)Gdx}, 



A(S?) 

for < to < £i < °° i we can write 

E[Y lm (K(ti))YL^ x At ))] (4.27) 

=y^(x)E[y, ro (»S(ti))] (4.28) 

=YC m {x)Y lm {x)E v (-m{h - t ) u ) by equation flUS]). (4.29) 

Thus, formula (|4.26p takes the form 

^QE[y m (^(ti))y^(^(t ))] =^Qy4(x)yi m (x)^(-^(t 1 -t n^(-/iK*i-*or). 

We also have that 

E[rU®*(*i))^m(»S(t2))] =E[y m (*Bf(t 1 ))y4(Q3f(t 2 ))], V<7G,SO(3) 
and this prove the first statement. 

We consider the case in which x ^ y: 

ii) From the assumption that the processes Q3J(ii) and 93^ (£2) are independent if and 
only if x ^ y we obtain that 

E[Y lm (K(h)) Y^Um = E[Y lm (K(h))] • E[y4(^(t 2 ))], iSx^y (4.30) 

where, from the formulae (|4.16[) and ()4.17|) of the Proposition 14.31 we have that 

E[y m (Q3£(ti))] = E v {-m (h - t Q f)Y lm {x) (4.31) 

and 

E[Y* m (mt 2 ))] = E v (-m (h - toTKMv). (4.32) 
For the sake of simplicity, we write 

= E v (-fii (ii - t ) u )E u (-ni (t 2 - t ) u ). 
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Prom ([PT]) and ([432"]) . the formula flg35D takes the form 

Im l>0 \m\<l 

By taking into account the addition formula (|3.6p we get that 

J2 Q E ^(^(t 2 ))] =£ C < £|(«o, ti, t2)-^i^((x, y» 

and the claimed result immediately follows. □ 
Remark 4.6. From (|4.ip roe arrive at the representation 

where the projection of%\ on Hi, space generated by {Yi m : |m < Z|}, written as 

T M( X ) = E OJm>!m(®J(t)), * > 0, Z > 
\m\<l 

is a random field representing the I th frequency component ofZ". The asymptotic properties 
ofZY t (as I — > oo) turn out to be very important in the analysis of the CMB radiation (see 
123[ \^4V- In particular, we have that 

|mi|<Zi |m,2|<i2 

From (|4.4p we get that 

E[^ ti (x)^ t2 (y)] =44 QE^^^))^^^))] 

|m|<i 
\m\<l 

where in the last step we used the fact that processes *Bf,(ti) and 23^ (£2) are independent if 
and only if x 7^ y. From Proposition \4-3\ we get 

mLtA^ktM] = s h s h E °i M-mih -t r)EA-Mt2-t r)Y lrn (x)Yr m ( y ) 

\m\<l 

and, from (|3.7|) we arrive at 

2/4-1 
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which is the (space/time) covariance of the I th frequency component of1%. Formula (|1.13D 
can be therefore rewritten as follows 

E[S£(s)S£(y)] =^E[^ tl (x)T^ 2 (y)]. (4.34) 

l>0 

We are now interested in studying the high-frequency behavior of such covariance. Usually 
the angular power spectrum of T is assumed to be as C\ = l~ a g(l) where a > 2 (to ensure 
summability) and g(-) is a smooth function which converges to a constant as the frequency 
I — > oo (for instance g(l) = G\(l)/G2(l) and G\, G2 are polynomials of the same order). 
The important fact is that E u has different behaviour depending on v as we know. Let 
us take to = for the sake of simplicity. The symbols S\ 5\ in (|4.33[> say that the Ith 
components are uncorrelated for different frequencies whereas, for large I and every 
spherical distance (x, y) 

n%M+hi x )%M (f )] ~ Ct^r°l h~ v , v G (0, 1) if t = 
where 6\ = a + 3 > 2 and, for u = 1, 

where Q\ = a ~ 1 > 1 an d C,D > are some constant that do not depend on I or h. This 
means that series (|4.34p converges with rate of convergence depending on v and therefore, 
the Ith frequency component (random filed) ZY t has different covariance structure as I 00. 
Also, for v G (0,1), the time covariance decays more slowly than an exponential decay (for 
v=\). 

We next give the proof of Theorem 11.81 

Proof of Theorem [TT8l For v = 1, we have a Markovian diffusion as the Chapman- 
Kolmogorov relation (|3.27p entails. Due to the Markov property we have that 

lP{*B x (t2) € dz 2 , ^(h) G d Zl } =P{<B x (i 2 ) G dz 2 \ *B x (ti) e dz x } 

x W{<B x (t{) G tfei|58*(to) G dx} 
x F{% x (to) G dx} 

for t 2 > t± > to > where, for all x £ Sf, 

F{& x (t ) G dx} = 1/A(S?). 
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Thus, we obtain that 

= exp(-fn(t 2 - h)) / Y lm (z 1 )Y l * m (zi)u 1 (zi,t 1 ;x,t )\(dzi) 



exp(-/i Z (t 2 - *i)) ^ exp(-/i 72 (ti - t )) 

72^2 

x / Y lm ( Zl ) YC m (zi)Y 12K2 (z{)Y* {x) \{dz x ) 
exp(-/^(i 2 - h)) exp(-/x 72 (ti - t )) Y* 2K2 (q 



72 «2 



(2Z + 1) 2 (2 72 + 1) ( M 72 W i / 72 



47r VOOOyVm — m k 2 

where we have used the formula (|A.10j) and the fact that Y^ m {z\) = (— l) m Y/_ m (zi). We 
write 



ilm _ I lV n /o; , i\ 2 ^2 + 1 (I I 72 \ I I I 72 



C™K 2 ~ ( I)™ (2^ + 1) y 

and therefore we get that 

ElY^^Y^ih))} = e-^-^J2 C ^2 e-^-^Y; 2K2 (x). 

72 «2 

Prom [3D] 

I I 



o o o ; 

and the fact that (see formula (|A.6j) or the book by [3]) 

+/ 



£(-W Z Z 72 )=^ K2 V2TTT 

^ Y m — m k 2 J 72 2 



E/^m = 2/ + 1 rOO 
72 K2 a/47t 72 K2 ' 

rn=—l 



m=—l 

we obtain that 



By recalling that Yqq(x) = 1/V47T we get that 

9/ 1 1 

^EPWCa-Cti))^^^))] = exp(- W (t2 -ti)) (4.35) 

m 

which does not depend on x and to- From this and formula (|4.26p we arrive at the claimed 
result. □ 
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Remark 4.7. The (equilibrium) covariance function of the TRD is given by 

4ir 1 
CoVgjfo u) := —E[Y 10 (<B*(i + h)) Y 10 (<8£(t )) = ^E v {-m h»), h> 0. (4.36) 

Indeed, from the fact that 

P{<B*(i) G dy,<B*(t ) G dx} =P{<B*(t ) G dx}P{Q3*(i) G dy|»J(t ) G dx} 

= j^n® X v(t) £ dy\<B*(to) £ dx}, 

we have to consider the integral 
1 



cos t? sini? cos i?o sini?o t P- l t; $0; ¥>0) *o) ^ di?o- 



5y taking into account the formula (see I40\l ) 

Y 10 {#,<p) = J-jLcaa# (4.37) 

V 47T 

and Proposition \4-4\ for I = 1 and m = 0, the above integral becomes the covariance (|4.36p . 
Cov« s (/ i ; !/) = ^E[Y 10 (®*(t + /i)) Y w (K(to)) ■ 



Also, for v = 1 we have that 



Cov s (h;l) = -exp(-/ii /i) 



which means that the rotational Brownian motion has an exponential memory kernel whereas, 
for v G (0, 1), formula (|4.36p says that the TRD has a covariance function with polynomial 
memory kernel, in particular Covr$(h; v) ~ h~~ v as h — )■ oo. Thus, we obtain that 

oo 

^Cov>b(/i;^) = oo (4.38) 
and therefore the TRD has long range dependence for v G (0, 1) whereas, 

oo 

£ctov»(M) = 3^:731 ( 43g ) 

implies short range dependence for the rotational Brownian motion (v = 1). 
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A Background on Clebsch-Gordan coefficients 



The Clebsch-Gordan (CG) coefficients are used in mathematics and in particular in rep- 
resentation theory of compact Lie groups. In physics, CG coefficients represent a set of 
numbers arising in angular momentum coupling under the laws of quantum mechanics. In 
this work we deal with CG and therefore Wigner coefficients because of their useful con- 
nection with integrals of spherical harmonics. For a deep discussion on spherical harmonics 
and Wigner coefficients we refer to the book [3D] while, for the interesting technique related 
to the graphical theory of angular momentum we refer to the book (JJ. Here we only recall 
some relation which turns out to be useful in the text. We write the CG coefficients by 
means of the Wigner 3j-symbols as follows 



W3 m 3 
hmil2m,2 



(-1) 



h-h+mz 



V^h + 1 



h 

771,1 



h 

m 2 



h 

-m 3 



h h h 

m\ 771,2 m 3 



"I) 



(h + h + h + l) 



(k + m 3 )l(l 3 -77i 3 ) 



1/2 



(h +mi)!(Z 1 - mi)\{l 2 + m 2 )\{h - m 2 )\ _ 

{-l) z {h + h+m 1 - z)\{h ~ mi + z)\ 



E 



z\(l 2 + h-h~ z)Kh +m 3 - z)\{h -h-m 3 + z)\ 



(A.l) 



with \rrii\ <U,i = 1,2,3. The Wigner's coefficient can be explicitly written as 

~(h +h- h)Kh -h + k)K-h + h + hV- 1 1/2 



where the summation runs over all z's such that the factorials are non-negative. For the 
3j-symbols we have the following properties. 

Triangle conditions: the Wigner's 3j coefficients are real-valued and they are different from 

zero only if mi + 7712 + 7713 = 0. 

Parity: 



^1 h h 
mi 7772 Tn 3 

Symmetry: for all h,hh > 



(-1) 



h+h+h 



h 

-mi 



h 

-777 2 



h 
-7773 



(A.2) 



h h h 

mi 7772 m 3 



O rt honor mality : 



E 



mim2 



h h h 

7772 m 3 mi 



h h L 

mi 7772 M 



h h a 

mi 777 2 W 



^3 ^1 h 
7773 777i m 2 



(2L + 1)' 



(A.3) 



(A.4) 



I I 



m 



-771 



7 

K 



5° 5° V2T+T, 



(A.5) 
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and 



where 



y^(2Z + 1) f h h 1 )( \\ \l 1 )=S£ffi (A.6) 

Ira 

y f h h h\f h h h \ = 1 (A7) 

^— ' \ mi 777,2 ?7i3 / V mi m2 W3 / 



mim2m3 



is the Kronecker's delta symbol. 



For the integrals involving spherical harmonics we recall that 



and 



/ \y l V d& BrntiY^fotfY^frtp) = (-ir 2 ^5"r (A.9) 
Jo Jo 
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dip d#sm0Y hmi (#,<p)Y hr n(0,<p)Y h ^(0,<p) (A.10) 
o Jo 



(2Zi + l)(2Z 2 + l)(2Z 3 + l) / h h h \ ( h h h 



47r \ / \ mi ni2 
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